function y=multvt(mu,C,v); 
global dim 
y=((randn(1,dim)*C)./(sqrt(gamrnd(v/2,2,1,dim))))'+mu;
% ===================================================
% function y=multvt(mu,C,v); 
% 
% Draw a sample y from a multivariate t-Student 
% t_p(mu, (C'*C)^(-1) , v ) 
% Following the notation in Lubrano et al 
% 
% Mean(y)= MU 
% Var(y) = (C'*C)/(v-2) = Sigma/(v-2) 
%          i.e. C'*C=Sigma  & C=chol(C) 
% See MULTRAND.m for a different parametrization following 
% Gelmand et al. 
% Alejandro Justiniano  December 20 2006 
% ==================================================
